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ABSTRACT 

We report the discovery of a broad emission feature at ~ 0.7 keV in the spectra of the ultra- 
compact X-ray binary 4U 1543-624, obtained with the high-resolution spectrographs of the 
XMM-Newton and Chandra satellites. We confirm the presence of the feature in the broad 
band MOS2 spectrum of the source. As suggested before in the literature, the donor star in 
this source is a CO or ONe white dwarf, which transfers oxygen-rich material to the accretor, 
conceivably a neutron star. The X-rays reprocessed in this oxygen-rich accretion disc could 
give a reflection spectrum with O VIII Lya as the most prominent emission line. Apart from 
the feature at ~ 0.7 keV we confirm the possible presence of a weak emission feature at ~ 6.6 
keV, which was reported in the literature for this data set. We interpret the feature at ~ 0.7 keV 
and ~ 6.6 keV as O VIII Lye and Fe K a emission respectively, caused by X-rays reflected 
off the accretion disc in the strong gravitational field close to the accretor. 
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1 INTRODUCTION 

The ultra-compact X-ray binaries (UCXB) are a subclass of low 
mass X-ray binaries (LMXB), which have very short orbital pe- 
riods, typically less than 80 min. In this type of systems the 
hydrogen-poor donor star transfers material to a black hole or a 
neutron star. The UCXB candidate 4U 1543- 624 has an orbital pe- 
riod of 18.2 min (Wang & Chakra bart\l2004h . The source is placed 
near the Galactic plane (/ = 322°, b = -6°) at a distance of ~ 7 kpc. 
It is a low inclination system (i < 70°), since it lacks eclipses or dips 
in the light curve. Type I X-ray bursts have never been reported for 
this source, therefore the nature of the compact object is uncertain. 
The donor star is most probably a CO or ONe white dwarf, based 
on the presence of ox ygen and carbon lines in the optical spectra 
dNelemans et alj2 006') and neutral neon in the X-ray spectra of this 
source with abundance wel l above the abun dance expected in the 
interstellar medium (ISM) jjuett et al.ll200lh . 
The X-ray photons, that are reprocessed in this oxygen-rich mate- 
rial transferred to the accretor may show up in the spectrum in the 
form of a broad emission O VIII Lyo- line (at ~ 0.65 keV), as re- 
cently found using X MM-Newton data in a similar oxygen-en riched 
UCXB 4U 0614+091 dMadei et alj20ld) . ISchulz et al.l ( l2010h con- 
firmed the presence of the OVIII Lyor line in 4U 0614+091 us- 
ing Chandra data. This line is the second strongest line after Fe 
Ka (at ~6.4 keV) in the reflection spectrum in the case of an ac- 
cretion disc with solar abundances. However, it may become the 
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most prominent line, when a significa nt overabundance of oxygen 
in the accretion disc is considered dBallantvne et al.ir2002t) . Since 
the line is thought to originate from the inner part of the accre- 
tion disc it is broadened by the effects of strong gravity close to 
the ac cretor, like gravitat ional redshift and the relativistic Doppler 
effect dFabian et alJ IT989). Additional broadening may come from 
Compton upscattering on hot electrons in the ionized disc material 
(Ball antvne et al.ll20021) . 

Evidence for an emission feature at ~ 0.7 keV i n 4U 1543-62 4 
was found befor e in ASCA and BeppoSAX data dSchultdl2003r) . 
lJuett et al.l fcOOll) used ASCA data to show that the feature could 
be fitted with a higher-than-solar neon abundance in the line of 
sight. The analysis of the high-res olution data of 4U 1543-62 4 con- 
firmed the overabundance of neon jjuett & Chakrabartvl2003f) . The 
authors reported however a 10-20% residuals around ~ 0.7 keV 
(18 A), which they ascribed to an calibration uncertainty. 
Emission at around 6.4 keV indicating the presence of the Fe Ka 
line was detecte d for this sou rc e in the high state sp e ctra o f RXTE 
and EXOSAT dSchultzl [2003h . lJuett & Chakrabartj] d2003h found 
also a Gaussian feature in the XMM-Newton data, discussed in this 
paper. Th e param e ters o f the Gaussian are comparable to the one 
found by ISchultzl d2003l) in RXTE data. The authors however re- 
ject the interpretation of this feature as due to Fe emission, be- 
cause of the large width and instability of the XMM fits. A nalysis 
of the same XMM-Newton data set done bv lNg et al] ({2010) reveals 
the presence of a much nar rower Fe Ka line than that reported by 
lJuett & Chakrabartvl d2003h . 

We reanalyze the high-resolution and broad band spectra of 4U 
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15 43-624 obtained by the XM M-Newton satellite. Since the work 
of lJuett & Chakrabartvl J2003h . the calibration of the RGS instru- 
ment improved sufficiently to show that the residual at ~ 0.7 keV is 
not a calibration issue, but a real emission feature. Furthermore, the 
feature is not present in the sources of similar flux and spectrum. 
In the spectrum obtained with the HETGS instrument on board of 
Chandra satellite we find a residual at ~ 0.7 keV similar to the one 
in the RGS spectrum. We attribute this feature to a broad O VIII 
Lycf line. Additionally we investigate the presence of the emission 
feature at ~ 6.4 keV in both data sets. 

2 OBSERVATIONS AND DATA REDUCTION 

XMM-Newton observed 4U 1543-624 for ~ 50 ksec on 2001 Febru- 
ary 4 starting at 13:17 UT. During the observation the two Re- 
flection Grating Spectrometers (RGS) and the European Photon 
Imaging Cameras MOS1, MOS2 and pn were collecting data. The 
MOS1 and pn cameras were operated in the Timing mode, while 
MOS2 was operated in the Full Frame (Imaging) mode. Because 
of the larger uncertainties in the spectral calibration for the Tim- 
ing mode with respect to the Imaging mode of the MOS cameras 
we do not analyze data obtained from MOS 1 camera. In the case 
of pn data, it was reported that pn camera in the Timing mode is 
not well calibrated in the soft part of the spectrum. Therefore in 
this analysis the part of the pn spectrum below 1.5 keV is excluded. 
We reduce the data using the XMM-Newton Data Analysis software 
SAS version 10.0. The observation is not contaminated by soft pro- 
ton flares. We extract the MOS2 observation with pixel pattern < 
12. The MOS2 observation suffers from pile-up, therefore we ex- 
clude events from within a circle with a radius of 15 arcsec centered 
on the source position. We check using the epat-plot tool that there 
is no evidence for pile-up in the selected annulus. The background 
spectrum for the MOS2 observation is extracted from CCD-3. The 
net source count rate is 9.37 ± 0.01 c/s. We extract the pn observa- 
tion with pixel pattern < 4 in rawx from 30 to 45. The background 
is extracted in rawx from 10 to 20. The net source count rate is 
92.41 ± 0.05 c/s. The data collected by RGS are reduced using the 
standard software pipeline which generates source and background 
spectra as well as response files. The net source count rate for RGS 1 
and RGS 2 are 4.04 ± 0.01 c/s and 5.75 ± 0.01 c/s, respectively. 
Chandra observed 4U 1543-624 for -30 ksec on 2000 September 
12 starting at 08: 16. During the observation the High Energy Trans- 
mission Grating Spectrometer (HETGS) together with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) was used. The HETGS 
includes two transmission gratings: the Medium Energy Grating 
(MEG) and High Energy Grating (HEG) covering the wavelength 
range 2.5-31 A and 1.2-15 A respectively. We use processed data 
available on the Transmission Grating Catalog and Archive (TG- 
Cat, http://tgcat.mit.edu/!. The net source count rate for MEG is 

10.8 ± 0.2 c/s and for HEG is 5.25 ± 0.0 2 c/s. 

We fit the data using the xspec package dAmau 3 1 19961) . Errors on 
the fit parameters reported throughout this paper correspond to a 
68% confidence level for each single fit parameter (A^- 2 = 1). In the 
case of the Chandra observation we use the Cash statistics, since 
most of the spectral bins in the region where the spectral edges oc- 
cur (see Sec. 3.2) contain less than 30 counts. 

3 ANALYSIS AND RESULTS 
3.1 Continuum model 

The model reported in the literature for the XMM-Newton and 
Chandra data sets, analyzed in this paper, consists of an absorbed 



powerlaw and black body djuett & Chakrabartvl [2003). Another 
model includes an absorbed black body plus a disk black body 
dNeetall koiO) and w as used to fit on ly EPIC-pn data, in the en- 
ergy range 0.7-10 keV. ISchultzl (2003) used the model comprised 
of black body and Comptonization components to fit RXTE, Bep- 
poSAX and ASCA data. The resulting high optical depths (r > 10 
in most cases) of the Comptonizing plasma mean that the spectra 
are close in shape to a cut-off powerlaw. Since the limited energy 
range (until 10 keV) of EPIC and HETGS data (until 7 keV) does 
not allow us to constrain the parameters of Comptonization model, 
we choose combinations of black body and powerlaw components. 
We simultaneously fit the RGS, MOS2 and pn data obtaining fits 
with similar^ 2 , value using the following models. Model A: a bro- 
ken powerlaw (bknpower model in XSPEC) together with a disc 
black body, absorbed by cold circumstellar/interstellar matter (TB- 
new). The TBnew model is also used in the Model B and Model 
C. Model B: an absorbed disc black body (diskbb in XSPEC) and 
a black body (bbodyrad). Model C: an absorbed cut-off powerlaw 
(cutoffpl) and a black body. The best-fit parameters are listed in the 
Table 1 . In the case of HETGS data, we obtain the best fit to MEG 
and HEG data (fitted simultaneously) using an absorbed powerlaw 
and a black body (model D). We add a multiplicative factor (con- 
stant) in front of all of the models in order to account for possible 
cross-calibration uncertainties of the instruments. For the reference 
abundances of ele ments in the neu tral absorber we use the proto- 
solar abundances of Lodden] J2003h . 

The RGS data are fitted in the energy range from 7 A to 30 A, 
MOS2 in the range 0.4 keV to 10 keV, pn in the range 1.5 keV to 
10 keV, MEG in the range from 2 A to 26 A and HEG in the range 
2 A to 13 A. The source flux in the energy range 0.5-10.0 keV for 
the XMM-Newton observation is 7.37 ± 0.01 x 10 -10 erg/cm 2 /s (cal- 
culated using MOS2 data and Model C) and for the Chandra ob- 
servation is 10.02 ± 0.08 x 10~ 10 erg/cm 2 /s (calculated using MEG 
data). 

3.2 Measuring abundances in the RGS & HETGS data 

We start the analysis of the region around O VIII Lyo- (at ~ 18 As; 
0.7 keV) from fitting absorption edges prominent in this part of the 
spectrum (see Fig. 2). We measure the depth of the following ab- 
sorption edges: the Ne K-edge (at 14.3 A), the Fe L-edge (at 17.3 A) 
and the O K-edge (at 23.05 A) by fitting the spectra locally. We use 
the following wavelength ranges for these fits: 21 - 26 A for the 
O K-edge, 16 - 18 A for the Fe L-edge and 12 - 16 A for the Ne 
K-edge. For these local fits the continuum model is fixed to the one 
found from the overall fit to the RGS, MOS2 and pn data. 
In order to determine the depths of the spectral edges of oxygen, 
iron and neon in MEG and HEG data we use the same wavelength 
ranges as used for the local fits to the RGS data. We fix the con- 
tinuum model parameters to those found from the overall fit to the 
MEG and HEG data. The results are included in the Table 1 . 

3.3 Emission features at ~ 0.7 keV and ~ 6.4 keV 

We refit the RGS, MOS2 and pn spectrum with the column den- 
sities (abundances and N H ) of oxygen, iron and neon fixed to the 
values measured from the local fits to the RGS data. We fix the Nh 
to the value obtained from fitting the overall spectrum excluding 
the region from 16 A to 21 A. The fits for all investigated con- 
tinuum models reveal structure around 0.7 keV (see Fig. 1). In the 
case of model B and C for the continuum there is also an emission 
feature present around 6.6 keV. We do not detect a significant emis- 
sion feature at 6.6 keV while using model A. It is conceivable that 
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Figure 1. The continuum models fitted to the XMM-Newton: RGS1 (red) + RGS2 (green)+MOS2 (blue)+pn (black) and Chandra MEG (red)+HEG (black) 
spectra of 4U1534-624. The panels show the unfolded spectra and the model components (dotted lines). Top left: XMM-Newton data set fitted with model A: 
TBnew*(bknpower+diskbb), Top right: XMM-Newton data fitted with model B: TBnew*(bbrad+bbrad+diskbb), Bottom left: XMM-Newton data fitted with 
model C: TBnew*(cutoffpl+bbrad). Bottow right: Chandra data fitted with model D: TBnew*(pl+bbrad). The residuals show a broad emission feature at ~ 0.7 
keV and weak feature at ~ 6.6 keV w hen model B and C is used. The absorption features at ~ 1.84 keV and ~ 2.28 keV in the pn data are probably caused by 
calibration uncertainties <Ng et al.l2O10l) . 




Figure 2. Left: the continuum model C overlaid on the XMM-Newton RGS spectrum of the 4U 1543-624. Right: the continuum model D overlaid on the 
Chandra MEG spectrum. Both spectra show prominent absorption edges: O K-edge at 23.05 A, Fe L-edge at 17.3 A, Ne K-edge at 14.3 A. The residuals show 
a broad emission feature at ~ 0.7 keV. 



this is caused by the fact that the energy of the powerlaw break in 
this case is close to 6.6 keV. 

First, we model the emission features with a Gaussian. The cen- 
troid of the line for all of the continuum models is around 0.69 keV 
with a « 0.06 keV. As for the feature in the high energy part of the 
spectrum, the centroid is around 6.7 keV with <x =s 0.4 keV and 6.6 
keV with <r « 0.6 keV for models B and C, respectively. 
The major improvement of the xl value comes from the inclusion 
of the line in the soft part of the spectrum (Ax 2 = -1555, -1854 
for model B and C, respectively, see Table 2). Further, minor im- 
provement of the xl value comes from inclusion of the line around 



6.6 keV (Ax 2 = -49,-150 for model B and C, respectively). 
We fit MEG and HEG data with the abundances of oxygen, iron 
and neon fixed to the values obtained from the local fits and N H ob- 
tained from fitting the overall spectrum excluding the region from 
16 A to 21 A. The fit reveals the presence of an emission feature 
around 18 A as well. Inclusion of the Gaussian line improves the fit 
to MEG and HEG data by AC tM , = -101. The center of the Gaus- 
sian line is at ss 0.7 keV and the a w 0.06 keV. We do not detect a 
significant emission feature around 6.6 keV. 
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Table 2. Parameters of the Laor profiles fitted to emission features found at ~ 0.7 keV and 6.6 keV in the XMM-Newton and Chandra data. In the case of 
XMM-Newton data three different continuum models are used (see Table 1 and Sec. 3.1). The values of abundances are fixed to the values found form local fit 
and the value of Nh is fixed to the one found from the fit with excluded region from 16 A to 21 A. The energies of the lines E and standard deviations tr are in 
keV, equivalent widths in eV. 
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E 2 


f2 


EqW 2 


X 2 v /d.o.f. 


4* 


A 


0.691 ± 0.002 


0.055 ± 0.002 


27.2 ±0.9 


1.89/2920 
1.42/2917 


-1415 












B 


0.689 ± 0.002 


0.058 ± 0.002 


31 ± 1 


1.93/2920 
1.40/2917 


-1555 


6.73 ± 0.06 


0.37 ± 0.05 


26 ± 11 


1.38/2914 


-49 


C 


0.691 ± 0.002 


0.064 ± 0.002 


35 ± 1 


2.05/2921 
1.41/2918 


-1854 


6.55 ± 0.05 


0.61 ± 0.05 


67 ± 19 


1.36/2915 


-150 


D 


0.70 ± 0.01 


0.057 ± 0.008 


30 ±3 


Csm/d.o.f./N B iN 
1316/913/918 
1215/910/918 


AC 
-101 













Table 1. Parameters from the best fit continuum models to XMM- 
Newton (Model A, B, C) and Chandra data (Model D). Model A: 
TBnew*(bknpower+diskbb), B: TBnew*(bbrad+bbrad+diskbb), C: TB- 
new*(cutoffpl+bbrad), D: TBnew*(powerlaw+bbrad) The errors on the 
abundances of oxygen, iron and neon are calculated from the local fits. In 
this fit the region from 16 A to 21 A was excluded. 



Parameter 


Model A 


Parameter 


Model B 


N H [10 2I cnT 2 ] 


2.93 ± 0.02 


Nh [10 21 cm~ 2 ] 


2.54 ± 0.03 


Aoi 


1.11 ±0.02 


Aoi 


1.31 ±0.03 


ANel 


4.4 ±0.1 


ANel 


5.3 ±0.1 


A F eI 


0.82 ± 0.09 


A Fe I 


1.0 ±0.1 




1.861 ±0.007 


Tbbl [keV] 


0.515 ±0.004 


r 2 


3.36 ±0.03 


7**2 [keV] 


1.524 ±0.004 


E [keV] 


6.69 ± 0.03 


T dbb [keV] 


0.320 ± 0.004 


Tdbb [keV] 


0.441 ±0.005 


Xlld.o.f. 


1.41/2762 


Xlld.o.f. 


1.42/2761 






Parameter 


Model C 


Parameter 


Model D 


N H [10 2I cnT 2 ] 


2.62 ± 0.03 


Nh [10 21 cm~ 2 ] 


3.37 ±0.1 


Aoi 


1.25 ±0.03 


Aoi 


1.0 ±0.2 


AncI 


5.1 ±0.1 


ANel 


4.3 ± 0.3 


A F eI 


1.0 ±0.1 


A F eI 


0.9 ± 0.3 


r 


0.99 ± 0.03 


r 


1.71 ±0.01 


E cu , [keV] 


1.14 ±0.02 


7*4 [keV] 


0.221 ±0.008 


Tbb [keV] 


1.584 ±0.005 


C sta t/d.o.f./N BIN 


1110/815/821 


Xlld.o.f. 


1.46/2762 







3.4 Laor profile 

Based on the fact that the accretion disc is overabundant with oxy- 
gen and affected by the strong gravitational field of the accretor, 
a possible interpretation of the feature at ~ 0.7 keV is a rela- 
tivistically broadened emission line of O VIII Lya. Therefore we 
fit a Laor profile to the residuals present in the RGS+MOS2 and 
HETGS data. We fix the wavelength of the line to 0.654 keV (the 
position of O VIII Lyo-). The outer radius of the disc is fixed to 
1000 GM/c 2 . The fit to RGS+MOS2 data gives an inclination of 
~ 87° and an inner radius of ~ 4 GM/c 2 for model B and C and 
~ l3GM/c 2 for model A. 

Fitting the residual present in the HETGS data with the Laor profile 
gives a reasonable value of the inclination of 54° (see Table 3). 
Additionally we fit the Laor profile to both features at ~ 0.7 and 
~ 6.6 keV for model B and C, coupling the inner radius, inclina- 



Table 3. Parameters of the Laor profiles fitted to emission features found at 
~ 0.7 keV and 6.6 keV in the XMM-Newton and Chandra data. In the case 
of XMM-Newton data three different continuum models are used (see text). 
The energy of the line marked in italics is fixed during the fitting. 

Laor 

Model £[keV] i [deg] r in [GM/c 2 ] xlld.o.f. 

A 0.654 70 ±1.5 13 -°!o8 1.46/2916 

B 0.654 86.9 ±0.2 4.1 ±0.4 1.44/2916 

6.3 ±0.1 87.0 ±0.2 4.4 ±0.4 1.43/2914* 

C 0.654 86.7 ±0.1 3.8 ± 0.2 1.45/2917 

6.38 ±0.09 86.8 ±0.1 3.6 ± 0.3 1.42/2915* 

slat 

ld.o.f./N BIN 

D 0.654 54 ±1 3.1 ±0.7 1216/909/918 

* In this fit the Laor profile used for modelling the feature at ~ 0.7 keV is 
also included. 




E [keV] 



Figure 3. The best-fitting continuum model consisting of an absorbed re- 
flection component (dashed line), cut-off powerlaw and black body (dotted 
lines) overlaid on the RGS and pn spectrum of the 4U 1543-624. The upper 
frame shows the spectrum and model components, the lower frame shows 
residuals in the units of the standard deviation. 

tion and emissivity index. The energy of the line in the soft part of 
the spectrum is fixed to the energy of OVIII Lyo-. The energy of 
the line at ~ 6.6 keV we keep free due to the possible contribution 
from a number of iron ions (Fe I-Fe XXVI). In the case of both 
models B and C we obtain a result consistent with the case of the 
fit of a single Laor profile to the 0.7 keV feature for the inclination 
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Table 4. Parameters from the best-fit model consisting of relativistically 
broadened reflection component, cut-off powerlaw and black body, fitted to 
the RGSl,2+MOS2+pn data. Parameters fixed during the fitting are marked 
in italics. 



TBnew*(kdblur*reflionx+bbody+cutoffpl) 


N h [10 21 cm~ 2 ] 


2.64 ± 0.04 


A i 


1.26 ±0.02 




5.4 ±0.1 


A F eI 


0.66 ± 0.09 


r 


0.50 ± 0.05 


Ecu, [keV] 


0.9 ±0.1 


T bb [keV] 


1.481 ±0.007 


qkdblur 


2.77 ± 0.03 


r„, [GM/c 2 ] 


6.0 


i [deg] 


54.3 ± 0.4 


Lteflionx 


1.97 ±0.04 


Fe 


0.5 


f [erg cm/s] 


250 ± 11 


xlld.o.f. 


1.44/2913 



* Abundance of iron in the refionx model. Abundances of other elements 
in this model are fixed at solar values. 

and inner radius. 

Similar to the case of 4U 0614+091, the lines, apart from relativis- 
tic broadening, can be additionall y affected by the C ompton scat- 
tering in the ionized disc material dMadei et al.l2010h . which, if not 
taken into account, could artificially increase the inclination in the 
Laor profile fit. 

3.5 Reflection model 

We use the RGS, MOS2 and pn data to fit a relativistically 
broadened reflection model reflionx, available in xspec, which 
takes into account effects of Compton scattering in the ion- 
ized reflecting material. Relativistic broadening is modeled us- 
ing kdblur and for the Galactic a bsorption we choose again 
the tbnew model dWilmsetal.ll2007b . For the continuum model 
we choose the cut-off powerlaw and the disc black body (TB- 
new*(kdblur*reflionx+bbody+cutoffpl), see Fig. 3). The only ad- 
justable abundance in the reflionx model is the abundance of iron. 
Therefore, in order to mimic the enhancement of oxygen we reduce 
the abundance of iron by a factor of two. During the fit abundances 
of interstellar oxygen, neon, iron and N H are kept free. We obtain 
a fit with ^j; = 1.44 for 2913 d.o.f. (see Table 4). The inclination is 
~ 54° and the ionization parameter is ~ 250, which indicates that 
the line is broadened by Compton scattering as well. 

4 DISCUSSION 

We report the discovery of a broad emission feature near 0.7 keV 
in the high-resolution RGS, broad band MOS2 and HETGS data 
of 4U 1543-624. A similar feature was found recently in the oxy- 
gen e nriched UCXB 4U 0614+091 dMadei et alJioTol i ischulz et al.l 
l20ld) . A reasonable interpretation of this feature is that of a rela- 
tivistically broadened O VIII Lya line, which would be the most 
prominent line in the reflection spectrum in the case of an oxygen- 
rich accretion disc. Fitting a Laor profile to XMM data gives a high 
inclination, inconsistent with the lack of eclipses and dips in the 
light curve of the source. The line is most probably also affected 
by Compton scattering in the ionized material of the disc. Fitting 
a relativistically broadened reflection model we obtain a more rea- 
sonable value of the inclination of ~ 54° and significant Compton 
broadening. 

We have checked the spectra for the presence of an emission feature 



at ~ 6.6 keV and we have found a Gaussian shaped feature, when 
the continuum is described by the model that includes a disk black 
body and two black body emitters or black body and cut-off pow- 
erlaw. The centroid of the feature is at ~ 6.6 keV with cr ~ 0.5 keV , 
which is consistent with the values obtained by iNg et alj ( T2O10h . 
The energy centroid indicates emission from iron ions up to Fe 
XXV. However, the value of the energy found from fitting a Laor 
profile is ~ 6.4 keV, which suggests emission from iron ions from 
Fe I to Fe XIX. We do not detect significant emission at the posi- 
tion of the Fe Ko- line when disk black body plus broken powerlaw 
fit function is used. The energy of the powerlaw break however 
falls in the range of the iron K line emission, potentially masking 
its presence. We do not detect any significant emission feature at 
~ 6.6 keV in HETGS data . The same conclusion was reached by 
lJuett & Chakrabartvl ( T2003T) . 

Concerning the ISM/CSM absorption in the line of sight , we find 
an ove rabundance of neon in the absorber with respect to Lodders 
d2003l) in RGS and HETGS data, which corroborates the pr evi- 
ous measurements djuett et al.ll200ll ; |juett & Cha krabartv 2003). A 
residual around the oxygen edge suggests that oxygen may be also 
in molecular form, for which the fine structure of the absorption 
edge can be signi ficantly different compared to the atomic oxygen 
dPinto et alj|201fj) . The residual near the neon edge may indicate 
the presence of war m ionized gas in the line of sight (similarly to 
the case of Cyg X- 1 Jjuett et alj2006f) . 
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